The possibility of absolute calibration of analog detectors by using parametric 

down-conversion: a systematical study 



Giorgio Brida 1 , Maria Chekhova 2 , Marco Genovese 1 , Alexander Penin 2 , Ivano Ruo-Berchera 1 
1 Istituto Elettrotecnico Nazionale Galileo Ferraris, Strada delle Cacce 91, 10135 Torino, Italy 
2 Physics Department, M.V. Lomonosov State University, 119992 Moscow, Russia. 

In this paper we systematically study the possibility of the absolute calibration of analog photo- 
detectors based on the properties of parametric amplifiers. Our results show that such a method 
can be effectively developed with interesting possible metrological applications. 

PACS numbers: 42.50.Dv, 42.62.Eh, 03.67.-a; 



in 
o 
o 

(N 

> 
O 

0\ 



> 
m 

On 
O 



in ■ 

o ; 
^ ■ 

Oh! 
> ■ 

C ■ 

s : 
cr 



X 



I. INTRODUCTION 

Accurate calibration of photodetectors both in ana- 
log and in photon-counting regime is fundamental for 
various scientific applications, which range from "tra- 
ditional" quantum optics pj to the studies on founda- 
tions of quantum mechanics |2|], quantum cryptography 
0, quantum computation 4], etc. 

In traditional optical radiometry primary standards 
are based on absolute sources or detectors Syn- 
chrotron and blackbody radiations are absolute sources. 
The spectral radiance of radiation emitted by electron 
storage rings can be predicted using the Schwinger equa- 
tion provided that the magnetic field and the electron 
current have been measured. In a similar way, the radiant 
flux emitted by a blackbody at a known thermodynamic 
temperature can be calculated from Planck's equation. 
Both sources provide radiation with continuous spectral 
distribution, the dominant contribution being in the soft 
X-ray and UV spectral region for synchrotron radiation 
and in the IR spectral region for blackbody radiation. 
The relative uncertainty of both these sources is about 
1 part in 10 3 . Among the absolute detectors, there exist 
the following two types: thermal detector, called elec- 
trical substitution radiometer (ESR), and semiconduc- 
tor photon detector. The ESR is based on the electrical 
substitution principle, that is, the heating effect of the 
unknown optical radiation is compared with the heat- 
ing effect produced by a measured quantity of electrical 
power (Joule effect). The operation of these detectors 
at cryogenic temperatures [6j, below 20 K, allows to re- 
duce measurement uncertainty down to 1 part in 10 4 . 
The operation principle of semiconductor photon detec- 
tors underpins on the photoelectric effect, that is, the 
generation of free electron-hole pairs at the absorption 
of a photon. The quantum efficiency is defined as the 
average number of free electron-hole pairs produced per 
incident photon. In a high-quality silicon photodiode the 
quantum efficiency for the visible range is close to unity 
to within a few tenths of one percent. The deviation of 
the quantum efficiency from unity could be determined, 
independently from other radiometric measurements, by 
self-calibration technique 7] . Uncertainties of a few parts 
in 10 4 appear to be the limit of this technique with com- 
mercial photodiodes. There is a perspective of improve- 



ment to 1 part in 10 6 or better with custom photodiodes 
operated at liquid nitrogen temperatures § 

For what concerns single-photon detectors, classical 
calibration techniques are based on the use of a strongly 
attenuated laser source whose (unattenuated) intensity 
has been measured by means of a power-meter. The un- 
certainty of this kind of measurements is principally lim- 
ited by the uncertainty in the calibration of the very low 
transmittance required for reaching single-photon level. 

This limitation has initiated the study of an alter- 
native scheme, based on the use of photons produced 
by means of spontaneous parametric down-conversion 
(SPDC) , where photons are emitted in pairs strongly cor- 
related in direction, wavelength and polarization. Fur- 
thermore, photons of the same pair are emitted within 
tens of femtoseconds. Since the observation of a photon 
of a pair on a certain direction (signal) implies the pres- 
ence of another one on the conjugated direction (idler), 
when this last is not observed this occurs because of the 
non-ideal quantum efficiency of the idler detector, which 
can be measured in this way 0, 0, fill [l2L Il3j . This 
absolute technique (and others related [lj]) is becoming 
attractive for national metrological institutes to establish 
absolute radiometric standards because it relies simply 
on the counting of events, involves a remarkably small 
number of measured quantities, and does not require any 
standards. 

Because of the success of the SPDC scheme for cali- 
brating single-photon detectors, it is important to anal- 
yse if similar absolute calibration methods could be de- 
veloped for analog ones, eventually allowing the develop- 
ment of a calibration scheme operating in both regimes. 

A seminal attempt in this sense was done in [l5| follow- 
ing the theoretical proposal of [ljj. Nevertheless, these 
results were limited to the case of very low intensity (as 
we will show in detail later) and were very far from being 
developed to the metrological level. An accurate analysis 
of the possibility to calibrate analog detectors by using 
SPDC overcoming these limits is therefore demanded. 

Incidentally, the quantum efficiency r\ of analog detec- 
tors appears also in equations describing suppression of 
photon noise in parametric down-conversion using the 
feedforward [22| or feedback [21j,[23] transformations. 
Hence, such experiments could eventually be used to de- 
velop an alternative scheme for analog detector calibra- 



2 



tion. 

The purpose of this paper is a systematical study of the 
possibility to calibrate analog detectors by using para- 
metric down conversion. 

In section II we will give a short presentation of the 
SPDC scheme used in the photon-counting regime. In 
Section III we will analyze multimode SPDC, following 
a theoretical description of SPDC |16( developed with 
an account for the high-gain regime (of which we give a 
short summary). After discussing some general results, 
we suggest possible calibration methods and show how 
the scheme of Refs. |0, ^| can be derived as the low- 
intensity limit of one of them. Finally, in Section IV we 
will consider the possibility of calibrating analog detec- 
tors using the schemes for photocurrent fluctuation sup- 
pression. In particular, two-mode squeezing and feedfor- 
ward techniques will be discussed. For these schemes, we 
will analyse new possibilities and limitations. 



II. THE SPDC SCHEME FOR CALIBRATING 
SINGLE-PHOTON DETECTORS 

The scheme for calibrating single-photon detectors by 
using SPDC is based on the specific properties of this 
process, where a photon of the pump beam (usually a 
laser beam) "decays" inside a non-linear crystal in two 
lower- frequency photons, 1 and 2 (conventionally dubbed 
"idler" and "signal"), such that energy and momentum 
are conserved (uj pump — loi+lo 2 , k pump = ki + k 2 ). More- 
over, the two photons are emitted within few femtosec- 
onds. In synthesis, the calibration procedure consists 01 
of placing a couple of photon-counting detectors down- 
stream to the nonlinear crystal, along the direction of 
propagation of correlated photon pairs for a selected pair 
of frequencies: the detection of an event by one of the 
two detectors guarantees with certainty, due to the SPDC 
properties, the presence of a photon with a fixed wave- 
length on the conjugated direction. If N is the total num- 
ber of photon pairs emitted from the crystal in a given 
time interval and (Ni), (N 2 ) and (N c ) are the mean num- 
bers of events recorded during the same time interval by 
the signal detector, the idler detector, and in coincidence, 
respectively, we have the following obvious relationships 



(Ni) = mN; (N 2 ) = m N, 



(1) 



where 771 and 772 are the detection efficiencies in the signal 
and idler arms. The number of coincidences is 



(N c ) = V1V2N, 



(2) 



due to the statistical independence of the two detectors. 
Then the detection efficiency can be found as 



This simple relation, slightly modified by taking into 
account the background subtraction and corrections for 
the acquisition dead-time, is the basis for the scheme for 
absolute calibration of single-photon detectors by means 
of SPDC, which reaches now measurement precision com- 
petitive with traditional methods fl3l ] . 



III. ANALOG DETECTION 



A. Basic formulas 

In order to study the possibility of absolute calibration 
of analog detectors we are interested in a model of SPDC 
working at any values of parametric gain. Firstly, the rea- 
son is the necessity to work with quite large intensities 
yielding continuous photocurrent. Secondly, we need to 
explore new possibilities to use the properties of SPDC 
for calibration without the usual low-gain limitation. A 
theory suited for these purposes is developed in and 
literature cited therein. For simplicity of description we 
will consider type-I SPDC in the degenerate case where 
the signal and idler frequencies are lj% — u)2 — u)p '^ ap . 
In the limit of monochromatic and plane-wave pump ap- 
proximation, only pairs of modes with opposite trans- 
verse wave vectors, q and — q, and with frequencies io\ — SI 
and <x>2 + SI are coupled as a consequence of energy and 
transverse momentum conservation. We can write the 
input-output transformation relating the field operator 
a(q, SI) at the input face of the nonlinear crystal to the 
field operator b(q, St) at the output face: 



b(q, SI) = u(q, Sl)a(q, SI) + v(q, Sl)a' l (-q, -SI). 



(4) 



The coefficients u and v are considered in |17| . For our 
analysis, we are interested not in the exact form of u and 
v, but rather in the properties 

Kq,ft)| 2 -Kq,tt)| 2 = l, 
u(q, Sl)v(-a, -SI) = u(-q, -Sl)v(q, SI), (5) 

(6) 

which guarantee the conservation of the free-field com- 
mutation relations [6(q, SI), 6^(q, SI)] = <5(q— q')5(0 — SI') 
and [6(q,Q),&(q,Q)] = 0. 

The far field observed in the focal plane of a thin lens 
of focal length / is obtained in jn| from the near field 
by means of the following transformation: 



B(x,t) = t-7 / dx'b(x.',t)e-^ 

A s/ JSa 



—l < ,t x-x 



(7) 



rfi = (N C )/(N 2 



(3) 



where A s = is the central free-space wavelength of the 
down-converted light and Sa is the transverse area of the 
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domain where PDC takes place. In practical situations 
it can be identified with the effective cross-section area 
of the pump beam. We stress that Eq. J2J) is correctly 
usable only for the calculation of normally-ordered corre- 
lation functions because it does not conserve the correct 
commutation relations. According to (@J, (JJJ), one can 
write the far field as 

•[u(x', H)aAx!, Si) + v(x', Sl)J(-^-x', -SI)], (8) 

(9) 

where 

u(x,Sl) = u(—x,Sl), 

As/ 

v(x,Q) = v(—x,Q). (10) 

The spatial variation scale of these coefficients is on the 
order of xq = yj ^%f, l c being the length of the crystal. 
One can interpret :ro as the transverse coherence length of 
SPDC in the focal plane of the lens. The variation scale of 
u and v in frequency, let us denote it SIq-, represents the 
typical bandwidth of SPDC in the temporal frequency 
domain and T co h = jj- is the coherence time. At the 
same time, 

is the diffraction pattern in the far-field plane due to the 
finite transverse size of the system. Also here we are not 
interested in the exact form of p(x), which is determined 
by the shape of Sa , but we observe that its typical size is 
Sdiff = (A s /) 2 /Sa and its amplitude g^— . Hereafter let 
us assume the size of Sdi // to be much smaller than the 
coherence length of PDC in the far field, i.e., Xdiff <C xq. 
Within this approximation the integrals can be evaluated 
considering p(x) as a delta function. Now we can calcu- 
late the mean value of the photon flux density operator 
J(x, t) = B< (x,t)B(x,t) in the detection plane. By us- 
ing Eqs. (jHJ and JI), and considering the vacuum as the 
input state, one obtains 

(7(x,i)) =— ?— f dSl\v(x,Sl)\ 2 . (12) 
Sdiff J 

The physical meaning of this quantity is the mean 
number of photons crossing the detection plane at point 
x at time t per unit area and time. The integral func- 
tion is usually referred to as the spectral gain of SPDC 
and its height represents the mean number of photons 
per coherence time. 



We are also interested in the second-order correlation 
function of the intensity fluctuations defined as 

<«(x,i)«(x',0) ee (7(x,t)7(x',0)-<7(x,t))<7(x',t')), 

(13) 

where (7(x, t)I(x', t')) is determined by the joint proba- 
bility of a photon reaching the detection plane at x' at 
time t' and the other one, at x at time t. It is convenient 
to distinguish between the two contributions, one due to 
the autocorrelation inside one beam (signal or idler), the 
other one due to the cross-correlation between the two 
beams: 

(6I(x, t)6I(x',t')) = Qn(x, t, x', t') + £i 2 (x, t, x', t'). 

(14) 

Here Qn (x, t, x', t') represents the autocorrelation contri- 
bution, when the distance between x and x' is less than 
x diff- On the contrary, Qn{x.,t, x',t') is different from 
zero only if x' ~ — x within Xdiff- Therefore, it is the 
term responsible for the well-known cross-correlation be- 
tween signal and idler beams. From Eqs. (JSJ and J2J, it 
is easy to show that 

£n(x,i,x',0 - (7(x,t))<5(x-x')<5(t-t') + 
+ b(x-x')| 2 J J ^L [e -i{n-n') { t'-t) 

|£(x,n)| 2 |«(x,n')l 2 L 

(15) 

C? 12 (x,;,x',OHp(x + x')| 2 / J ^e-^n-nw-t) 

w*(x, n)w*(-x, -n)u(x, Sl')u(-x, -SI'). 

(16) 

The first term on the right-hand side of Eq. I|15fl is due 
to the commutator [B(x, t), B^(x', t')] = S(x—x')S(t—t') 
and it is responsible for the rise of shot noise during de- 
tection, while the second term is the normally ordered 
auto-correlation. In Eq. (|16fl the shot-noise term is not 
present because the commutator is null, being the cross- 
correlation calculated for x ~ — x'. Let us consider two 
detectors, 1 and 2, having quantum efficiencies ryi and 7/2, 
to register photons crossing two arbitrary but symmet- 
rically placed regions R\ and R2 in the detection plane. 
Let at least one of the two areas be much bigger than 
the diffraction area Sdiff [2(|- In fact, from the works 
on photon-counting detectors calibration it is well known 
that the area of the detector under investigation should 
be large enough to collect all the photons correlated with 
those incident on the trigger detector surface, otherwise 
the quantum efficiency is underestimated. Integration of 
Eq. l|12fl over i?i gives the photon flux reaching detector 
1: 

(70 - -L- f dx [ dSl\v(x,Sl)\ 2 (17) 

Sdiff jRi J 
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The autocorrelation and cross-correlation functions of 
the photon fluxes can be obtained by integrating defini- 
tion l|18[l. respectively, over R\ x i?i and R\ x i? 2 and 
using Eqs. (JTSJl and JT^: 

(7 1 (i)/ 1 (0) = (/i) 2 + (/i)^-0 + 
+_!_ [ dx [ [ ^MY [e -i(n-n' W - t) 

b drff J Rl J J 27T 

\v^,n)\ 2 \v(^n')\% 

(18) 



(Ii(t)h(t')) = {h)(h) + / dx / / 

e -i(n-n')(f-*)y ( Xj ft)S*(_ x , -fi)v(x, n')5(-x, -fi'). 

(19) 

In the quantum theory approach the shot noise arises 
from the commutation relations of the free-field opera- 
tors, see, for example, |19|. Any losses in the channels 
must be theoretically described, so that the correct com- 
mutators are preserved until the moment of the mea- 
surement, i.e., the photon-electron conversion inside the 
detector. Thus, to take into account the losses due to the 
non-ideal quantum efficiency, we consider any real detec- 
tor as an ideal one preceded by a beam splitter with the 
transmission coefficient equal to the quantum efficiency rj 
of the real detector. We substitute the field B(x, t) enter- 
ing the beam splitter with the transmitted field C(x, t) 
defined as 



C(x, t) = ^B(x, t) + v/I - ^V(x, t), (20) 

where V(x, t) is the field operator for the second input 
port of the beam splitter, which is assumed to be in the 
vacuum state. From this commutator-preserving trans- 
formation, it easily turns out how to regard sub-unity 
quantum efficiency: it is sufficient to replace, in the nor- 
mally ordered products, B(x,t) with y/rjB(x.,t) 0]. 

B. Feasibility of analog detectors calibration by 
measuring SPDC correlations 

For our calculations, we will consider a few millimeters 
non-linear crystal pumped by a CW laser. If the waist of 
the pump, identifiable with the transverse cross-section 
of the system Sa, is relatively large, namely on the order 
of millimeter or more, the real system fits the model of 
SPDC discussed in the previous section. 

In any detectors the absorbtion of a single photon, 
i.e., any detection event, generates an electric pulse in 
the current or in the voltage having a profile fit) and a 
random area q n at a random time of occurrence t n . In 
the analog process, with a large incident photon flux, we 



cannot distinguish between two different pulses because 
they overlap. Rather, the information about the statis- 
tics of light is carried by the continuous fluctuations of 
the photocurrent. We express the current (or voltage) as 
a superposition of many pulses [2(j : 

n 

In an ideal instantaneous-response photocell, all values 
q n are equal to the charge e of a single electron and 
f(t) ~ S(t). For real detectors, the time constant r p 
is finite, and typically r p ~ Ins or more. Also, when 
a photodetector is operated in the avalanche multiplica- 
tion mode, this process gives rise to an internal current 
gain. The statistical nature of the multiplication process 
gives an additional contribution to the current fluctua- 
tions m. 

Since the probability density of observing a photon at 
time t is related to the quantum mean value of the photon 
flux (1(f)), we calculate the average current as 

(h(t))=J2(<Hnf(t-t ln )) = 
n 

J dtMfit-hJihit)) . (21) 

Analogously, the quantum-mechanical second-order in- 
tensity correlation function is determined by the proba- 
bility density to have a photon detected at time t and 
another one at time t' , see Eqs. (|13ll . Therefore the 
correlation functions for the current fluctuations can be 
expressed as 

(h(t)ii(t + t)) = ^^qinQlmfjt - h n )f(t - t lm + t)) = 

71, m 

dhdt'Milfit - h)f(t - t[ + r)(/i(ii)Ii(ii)), 

(22) 

(ii(t)i 2 (t + t)) = y^(gi w g2ro/ft - h n )f(t - t 2m + t)) = 

n,m 

dtidt 2 (giga)/(* - h)f(t -t 2 + r)(7i(ti)/ 2 (i 2 )), 

(23) 

where the first equation is the electron current autocor- 
relation function for each detector, one registering the 
intensity of the signal beam and the other registering the 
intensity of the idler beam, while the second one is the 
cross-correlation function between the electron currents 
produced by the two different detectors. 

According to Eqs. (|21|l and (|17ll . the mean value of 
the electron current (the analog of Eq. in photon- 
counting regime) is equal to 

(il)=Vl{qi)(h)- (24) 
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Here the quantum efficiency has been taken into account 
according to formula l|20l) . with the substitution (Ji) — > 
r/i(/i). The factor (qi) is the average charge produced in 
a detection event. 

The current correlation functions can be obtained sub- 
stituting Eqs. Ijl 8JI and ljl 9JI . respectively, in (|22|l and 
(|23|l . We note that functions in the integrals of (|18|) and 
(|19|l have significant values for |Q — Sl'| ~ Slo- Roughly, 
this means that the correlation functions have a sinc-like 
behavior in time, with the central peak width equal to 
the coherence time of PDC T co h = n- ~ 10~ 13 s. As men- 
tioned above, the resolving time of a real analog detector 
is finite, and in general can be considered much larger 
than the SPDC coherence time. Thus any fluctuations 
in the intensity of light are averaged over t p during the 
detection process. So in the limit t p 3> t co /j we have 



(h(t)h(t + T)) = {t 1 ) 2 + m {qj)T(r) 
1 



(h) + Vi 



S d 



iff JRi 



,/x / |^(x,tt)| 4 



(25) 



and (the analog of Eq. (|3J in the photon-counting 
regime) 



(h{t)i 2 (t + r)) = (iiXia) + mV2(qi)(q2)H< 

bdiff jR t J 27T 

where we have defined 



(26) 



^(r)= / dtf(t)f(t + r). 



(27) 



The last two equations are the fundamental tools for 
studying the problem of absolute calibration of analog de- 
tectors and thus we are going to discuss them in detail. 
Despite (12511 and (|26|l seem to be quite symmetric, we 
observe some important differences. The presence of (ii) 
in the autocorrelation function is due to the shot noise 
contribution and for this reason the quantum efficiency 
r\ enters linearly, while in the current cross-correlation 
function the corresponding term is due to the high quan- 
tum correlation between the signal and idler beams of 
PDC and the quantum efficiency appears quadratically. 
It is equivalent to the right-hand side of Eq. (J2J for the 
counting regime and its presence is the key for absolute 
calibration. The second term inside the brackets, both for 
auto- and cross-correlation functions, is important only 
when the number of photons per coherence time is not 
close to zero and so the presence of two or more photons 
within that time is not negligible. In fact, the integral 
term can be estimated as v 2 (Ii) and can be neglected 
as long as v 2 <C 1, i.e., the mean number of photons 
per T co h is much smaller than one. Anyway, if the du- 
ration of the photocurrent pulse is much larger than the 
coherence time, this assumption does not prevent pho- 
todetection to be in a strongly analog regime, because 



a lot of photons can be absorbed during the pulse du- 
ration as well. The term proportional to (I) 2 is due to 
the presence of more than one photon in time r p and 
for that reason is more delicate. As we can observe in 
Eq. (|2*6*|> it can be neglected only if (I) <C F{t). Since 
the pulse f(t) has a height around l/r p , by Eq. I|27|l . 
max [Fir)] — J-(0) ~ 1/t p . So the condition becomes 
(I)t p <C 1, i.e., the number of incident photons during the 
resolving time of the detector should be much less than 
one, i.e. one should work in a non-overlapping regime. 
In principle, in this case one could distinguish between 
different pulses of the current and work in the counting 
mode. 

The usual definition of the quantum efficiency is the 
ratio between the number of photons detected and the 
number of photons incident on the detector surface. This 
definition is completely suitable in the case of counting 
detectors and is exactly the meaning of 77 in our paper. 
But in the case of analog detection, we cannot in prin- 
ciple distinguish between different current pulses. Thus, 
according to formula (|24|) , we adopt the definition of ana- 
log quantum efficiency as T = r)(q) — (i) / (/), having the 
meaning of the electron charge produced per single inci- 
dent photon, or the ratio between the electron flux and 
the photon flux. If the charge q produced per photon 
fluctuates it increases the current fluctuations. This ex- 
plains why (q\) appears in Eq. H25fl instead of (qi) 2 . In 
principle, the most general way to obtain an estimation 
of r\ working with the PDC light intensity in the photon- 
counting regime is dividing the coincidence counting rate 
(proportional to the cross-correlation function) by the de- 
tector counting rate (proportional to the intensity) . This 
method works because in the photon-counting regime the 
temporal shape of the current pulses and their width is 
not important; instead, one registers a single pulse or 
not registers it. This is not the case for analog detec- 
tion in which the pulse shape f{t) appears in formulas 
through the factor !F(t). In general, we do not know 
this function, and this makes the absolute calibration of 
analog detectors more difficult. However, as we are going 
to show, under some condition it is possible to overcome 
this drawback. 

Let us distinguish between three different regimes: 
very low intensity (I) , middle intensity (II) , high intensity 
(III). 

(I) (I)t p <C 1 (i.e., photocurrent pulses do not over- 
lap). For example, for a detector with a time constant 
t p = 10 ns, the corresponding photon flux must be below 
10 8 photons/s. In terms of power, for the wavelength of 
500 nm, it means about 10 pW. 

Eqs. (|2*5)l and l|2*6l become then 



(i 1 (t)i 1 (t + T))=r ll (q 2 1 )F(T)(I 1 ), 



(28) 



(ii(t)i 2 (t + T))=r h r ]2 {q l ){q 2 )F(T){h). (29) 
The same equations has been found in IQj and the 
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quantum efficiency has been estimated as 



T2 = ?72<<32> = 



(gl) / , w 

(«l) aWi; <*l(t)*l(* + T))- 



(30) 



This formula is not satisfying from the metrological point 
of view because of the presence of some unknown param- 
eter related to the statistics of charge fluctuations that 
we discussed previously and that has to be estimated 
in some other way. We suggest to avoid the problem 
by integrating Eq. (|29|) over time r. It could be done 
after the acquisition of the profile of the function has 
been performed. By definition 1|27|> . it is evident that 
J dr!F(T) = 1; integrating Eq. I^H} in r and dividing it 
by Eq. i|21fl. we obtain 



J dT{n{t)i 2 {t + T)) 



(31) 



As pointed out, another drawback of l|3U|l is the neces- 
sity to work only at very low intensity, where no overlap- 
ping between pulses happens. This, in terms of experi- 
ment, means that one has to work in the so-called "charge 
accumulation mode" , in which the electron charge is ac- 
cumulated until reaching some detectable threshold. In 
the case of avalanche devices, it corresponds to the pos- 
sibility of direct photon-counting. 

(II) (I)r p > 1 but still v 2 <C 1 (i.e., photocurrent 
pulses overlap but the parametric gain and photon flux 
are still quite low). Considering the same parameters as 
used in case I, coherence time r co h of the order of 100 
fs, and the requirement that v 2 < 0.001, this means a 
photon flux between 10 8 photons/s and 10 1 photons/s 
or power between 10 pW and 100 nW. 

Eqs. and |J2BJ| become 



(h(t)i 1 (t + T)) = (i 1 ) 2 +r h (q 2 1 )F(T)(h), 



(32) 



(ii{t)i 2 {t + T)) = (ii)(i2)+mV2(qi)(q2)F(T)(h). (33) 

We stress that in ^(J > to the best of our knowledge the 
only theoretical paper treating the absolute calibration 
of analog detectors using PDC, the very low parametric 
gain was assumed from the very beginning. In partic- 
ular, all terms proportional to v 4 were neglected. This 
corresponds to the approximation (I) in our treatment. 
Accordingly, in ^1 the limitations of Eq. (|30|l were not 
discussed. Our analysis shows that when the intensity is 
high enough to yield a strongly analog current, Eqs. 
and (|33|l should be used. Now we define the correlation 
functions of the current fluctuations as 

(6i k (t)6n(t + t)) = (i k (t)ii(t + T)) - 

- (i k (t))(ii(t + T)) (M = M). (34) 

We underline that these functions remain in principle 
experimentally estimable. So a new formula, similar to 



(|30ll . is available for analog quantum efficiency estima- 
tion: 



772(92) 



(q 2 ) {5h(t)5i 2 (t- 



0> 



{5n{t)5n{t + T))' 



(35) 



Once again, the drawback of this formula is the pres- 
ence of the unknown parameter (q 2 )/(qi) 2 that requires 
additional measurements to be performed. As before, 
integrating in r the expression for the cross-correlation, 
i.e., the definition Ij34(l for k = 1, 1 = 2, we obtain 



J dT{5i 1 {t)5i 2 {t + T)) 



(36) 



This equation represents one of the main results of the 
paper, since it shows that the absolute calibration of ana- 
log detectors by using SPDC is indeed possible. 

A drawback of quantum efficiency measurement in this 
regime could derive from the fact that the terms propor- 
tional to the square of intensity in (|32|l and (|33|l become 
rather large, much larger than the term proportional to 
the intensity - the one that provides calibration. We 
need to subtract this background, as it is done in 13411 . 
Although it is an easily estimable quantity, a little rela- 
tive uncertainty could generate a large uncertainty in the 
efficiency estimation: a limitative result when looking to- 
ward metrological applications. The physical reason for 
this behaviour can be found recalling that the quantum 
correlation so attractive in PDC has the scale of r co ^. In 
the analog regime this correlation is almost deleted be- 
cause of the averaging over a time r p ^> r co /j. Anyway, 
as long as we take (I)t p ~ 1 -j- 10, such a problem does 
not arise. 

An alternative could be working in the pulsed regime, 
in which the duration of any pump pulses is not so far 
from the coherence time of PDC and the distance be- 
tween them is larger than r p . Of course, to have a 
large number of photons during t p (strongly overlap- 
ping regime) we need to increase the parametric gain. 
A detailed study of this possibility will be presented else- 
where. 

(Ill) v 2 > 1 (i.e., high- intensity regime). 

In this regime each term of (|25|l and (|26[) is important 
and no general way can be found for the absolute cali- 
bration of analog detectors, at least with a CW pump. 
It can be shown that in the single-mode case the integral 
terms in (|25fl and (12(:i[l become equal to the first ones, 
proportional to the square intensity, and thus could be 
easily estimated. In the case of CW pump, single-mode 
detection requires very narrow filters and fast detectors, 
beyond realistic present technological possibilities. Fur- 
thermore, calculation shows that for analog calibration 
we need to know the transmission spectrum of the filters. 
Anyway, also in this case one can consider the possibility 
to reach single-mode detection by using a pulsed pump. 
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IV. SQUEEZING FOR CALIBRATION 



In the 1980's and the 1990's a lot of work has been 
done to demonstrate both theoretically and experimen- 
tally the possibility to obtain sub-shot-noise photocur- 
rent statistics taking advantage of strong PDC quantum 
correlation. The shot-noise level (SNL) is defined to be 
the lower limit to the photocurrent noise level, which is 
achieved for coherent states of the field. Basically, for 
what concerns two-mode squeezing, two different kinds 
of schemes have been used. The first one consists of de- 
tecting the currents from the two light beams (signal and 
idler) and subtracting them [24( . The variance of the dif- 
ference current or the difference between the numbers 
of generated electrons can be less than the same quan- 
tity measured for coherent beams. In the other scheme 
the information about fluctuations in one beam is used to 
manipulate the intensity of the second beam (feedforward 
technique) 22] , or directly the pump intensity (feedback 
technique) [23| in order to correct the fluctuations. The 
goal in this case is getting a reduction of photocurrent 
fluctuations below the shot-noise level for the detector 
measuring one of the beams. 

It is important that in both schemes the minimum 
reachable squeezing factor depends strongly on the quan- 
tum efficiency of the detectors. Thus, it is reasonable to 
consider the possibility to extract the quantum efficiency 
from the degree of squeezing. 

Let us first consider the two-mode squeezing introduc- 
ing the difference between the currents z_ = i\— 12. Using 
equations (|^3f> and it2fi|) the autocorrelation function can 
be easily evaluated. When we consider two balanced de- 
tectors, Tji(qx) = 772(92) = ^(9)1 collecting photons from 
symmetric and equal detection areas, and no fluctuations 
in the charge produced per photon occur, i.e., (q 2 ) = (q) 2 , 
we have 



(5i-(t)6i-{t + T)) = SNL(1 - r)), 



(37) 



where the shot noise level is given by SNL = 
2(q) 2 r)(I)J 7 (T). Let us discuss this formula. First of all, 
its validity does not disappear in the high-gain regime 
and at high light intensity because the integral terms 
of (|25|) and (|26|) cancel each other in the calculation of 
(|37[) . Quite surprisingly, it shows that some aspects of the 
quantum correlation, such as the possibility of shot-noise 
suppression, can be preserved when working with macro- 
scopic intensities, where the cross-correlation is domi- 
nated by classical terms proportional to the square of 
the intensity. Furthermore, it opens the possibility of cal- 
ibrating analog detectors with SPDC at high parametric 
gain. In particular, it can be useful in the calibration of 
detectors in which the electronic noise of the external am- 
plifier is dominant compared to the shot noise when work- 
ing at low intensities (cases (I) and (II) of the previous 
section). In this kind of detectors, like, for instance, sim- 
ple photodiodes, no avalanche or multiplication occurs 



and in general we can consider the condition (g 2 ) = (g) 2 
satisfied. So the requirements for the derivation of i|57|) 
are realistic in this case. It can be used for calibration di- 
rectly, by calibrating the SNL with the help of a coherent 
source producing the same average current as measured 
for one of the SPDC beams, or by integrating in r and 
then normalizing by (ii): 



J dT(Si-(t)5i-(t + T)) 



*i 



%)(!-??)• 



(38) 



From the experimental viewpoint, the validity of i|37|) 
is guaranteed if the detection volumes of the two detec- 
tors are conjugated. In other words, detector 2 has to 
collect exactly all and only the modes conjugated with 
those collected by detector 1. For that pertaining the 
spatial modes, detailed realistic calculation performed in 
as well as recent experiments show that the de- 
tection areas should be large enough to collect several 
spatial modes. By satisfying this condition, one would 
also increase the detected radiation power and hence, the 
signal-to-noise ratio. 

Finally, we note that an equivalent formula is valid 
for the counting regime without any assumptions except 

m = m = v- 



((SN-) 2 ) 
(N) 



2(1-7?), 



(39) 



where 7V_ = Ni — N2 is the difference between the num- 
bers of counts in detectors 1 and 2, and ((<5iV_) 2 ) is its 
variance. This method of calibration could be interest- 
ing because, unlike the method of coincidence counting 
based on Eq. (0), it only requires counting photodetec- 
tion pulses. 

Of course the assumption of balanced quantum effi- 
ciencies used for deriving l|37|) somehow complicates the 
aim of reaching the accuracy needed for metrological ap- 
plications, although the technique of balancing detectors 
is largely used in quantum optics. 

Finally, we consider the possibility to exploit the feed- 
back or feedforward schemes. In Ref [22|, the latter is 
studied in detail for optical parametric oscillator (OPO) 
above threshold where one has bright average intensity 
component and small fluctuations. It is shown that the 
minimum of fluctuations achievable can be expressed in 
terms of the frequency spectra of the fluctuations in each 
beam and the correlation spectrum in the absence of the 
feedforward action. Unfortunately the results of [23] can 
not be applied directly to SPDC because in the case 
of SPDC the resulting state is squeezed vacuum where 
high fluctuations of intensity occur. Anyway, descrip- 
tion of the feedforward scheme presented in |22fl can be 
applied to the currents. For instance let the current in 
the detector registering beam 1, ii, be varied by using a 
modulation mechanism conditioned by the measurement 
performed on beam 2. If i\ is the current registered af- 
ter modulation, in our time representation (for simplicity 
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with t = 0) the result of Ref. |22 

((^) 2 ) = ((^i) 2 > - 



| yields 

(5ii5i 2 ) 2 
<(^ 2 ) 2 > 



(40) 



This formula shows that is possible to measure the reduc- 
tion in the current fluctuations of detector 2 after modu- 
lation instead of measuring the cross-correlation function 
(5i 1 5i 2 )- Substituting (jS2Jl and into (gOl, we have 



- 1 - ?7lT72 / 2 \/ 2\ 



(41) 



If the two detectors have balanced quantum efficiencies 
and excess charge noise is absent, the above formula can 
be used for the efficiency estimation, representing a fur- 
ther interesting option in this sense. 

It is interesting to observe that in [2j| , the authors ob- 
tained the same suppression of fluctuations as given by 
(|41|l using a theoretical model describing a feedback pro- 
cedure applied to SPDC. The only difference is that our 
equation takes into account the statistics of the charge, 
which enters through the factor ,L,l{ . Even if the 

quantum efficiency is ideal, this factor prevents one from 
reaching perfect noise reduction. 



of performing calibration in the high-intensity regime of 
SPDC, where other methods fail, but with the limita- 
tion that one should have two detectors with balanced 
quantum efficiencies. Also, a possibility of estimating 
the quantum efficiency for the feedforward or feedback 
scheme of fluctuations suppression is considered, reach- 
ing analogous conclusions. 
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V. CONCLUSION 



Motivated by the necessity of a general absolute cal- 
ibration scheme for analog detectors for various appli- 
cations and aiming to extend the absolute calibration 
scheme from the single-photon to the analog regime, we 
have performed a systematical study on the possibility 
of applying SPDC calibration methods to the analog 
regime. Possibilities and limitations following from the 
specific properties of SPDC have been investigated. 

Our results show that measurement of the correlations 
in the output currents indeed can be used to extend the 
absolute calibration method to the analog regime, al- 
though the experimental implementation will require an 
accurate study and solution of some technical problems. 

In particular, it is shown that integration of the pho- 
tocurrent correlation functions in time allows one to 
avoid the measurement of the photocurrent pulse shape 
and to eliminate the necessity to know the statistics of 
electrons in the photocurrent. 

Also, our analysis showed the possibility to go beyond 
the regime of non-overlapping photocurrent pulses, which 
was used in earlier works, and to pass to higher intensi- 
ties. A possible way to reduce the background in the 
measured photocurrent correlation function, which will 
unavoidably accompany the transition to higher inten- 
sities, namely, passing to the pulsed regime of PDC, is 
outlined. 

Finally, we studied the new subject of squeezing as a 
tool for absolute photodetectors calibration. In partic- 
ular, two-mode squeezing is shown to be the only way 
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